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Abstract. We consider a variational problem in a planar convex domain, 
motivated by statistical mechanics of crystal growth in a saturated solution. 
The minimizers are constructed explicitly and are completely characterized. 



1. Introduction 

In understanding the physical phenomenon of droplet condensation or crystal 
growth, the central issue is to explain how a particular macroscopic shape of 
the growing droplet or crystal is determined by microscopic interactions of its 
constituent paricles. 

According to Gibbs' formulation of statistical mechanics, the probability of 
a microscopic configuration a is proportional to exp{—/3H{a)), where /3 > is 
the inverse temperature and H{-) is the Hamiltonian defining the energy of the 
system. Therefore the most probable configurations are the ones with minimal 
energy. In the "thermodynamical" limit of a large number of particles, this 
minimum becomes very sharp: the overall configuration of the system settles, 
up to minute fluctuations, to a well-defined deterministic structure. 

It turns out that the microscopic laws of atomic interactions give rise to a 
certain macroscopic quantity, the surface tension, which determines the droplet 
shape via minimization of the surface energy. Phenomenology of surface tension 
has been proposed by Gibbs in the late 1870's. In an important contribution, 
G. Wulff suggested in 1900 that for a growing crystal, its equilibrium shape is 
that of a ball in a metric generated by the surface tension (the Wulff shape). 

It has been furthermore observed experimentally that flat facets of a growing 
crystal may carry macroscopic but monomolecular "islands", whose shape is 
also determined by the surface tension. A mathematical approach to explana- 
tion of this phenomenon has been developed by S. Shlosman and collaborators 
in a series of works 114^181 [9], building upon his earlier work with R. L. Dobrushin 
and R. Kotecky [6]. 

A typical setting in this approach is represented by the following discrete 
model of crystal growth, which is a variant of the Ising model: fix an open 
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domain C of unit area and consider the three-dimensional lattice obtained 
by intersecting the cylinder O x [—1, 1] C with -^{^^ + (0, 0, |)), where N is 
a large integer parameter. At each node t of this lattice there is a variable at 
(the spin) taking values +1 (interpreted as "t belongs to the free phase") and 
— 1 (interpreted as "t belongs to the condensed phase"). The collection a = (at) 
is called the microscopic configuration of the system. 

Fix now the Ising Hamiltonian H{a) = — X^^j. "^s*^*' which describes 

a "ferromagnetic" interaction between nearest neighbors (equal values have 
smaller energy than opposite ones), and consider the canonical probability dis- 
tribution p{cr) = exp{—(3H{a))/Z. Here the normalization coefficient Z = 
exp(— /3i?(c7)) is defined by summation over all configurations that satisfy 
the so-called Dobrushin boundary condition: spins at outermost nodes {x,y,z) 
of the lattice have values if z > and —1 if z < 0. 

It turns out that in the limit of large N the main contribution to probability 
comes from configurations where the lower and upper halves of the lattice are 
filled respectively with — I's and -|-l's. In this equilibrium state, the numbers 
of -|-l's and —I's are asymptotically equal, so that Sn = Ylt^t ~ 
fluctuations of the flat surface dividing the two phases are logarithmic in N. 

A more interesting situation occurs when, in addition to the Dobrushin 
boundary values, the system is conditioned to have macroscopically more —I's 
than -|-l's: 

Sn = o"t = —mN"^ 
t 

with m > 0. In this case, depending on the value of m, the most probable 
state of the system may feature one or more monomolecular layers on top of 
the surface z = in the box x [— 1, 1]. A detailed account of the observed 
equilibrium states as m changes can be found in j9]. 

As proved in [13] , the behavior of this model in the continuous limit N ^ oo 
is closely related to the following variational problem: given an open set C 
and a value m G [0, +oo), find 
(1.1) 

mini / ip*{Du): ueBV{W), n = on M" \ 0, n(-) G N, / udx = m\, 
[Jn Jn } 

where (p* is some given general norm on M". Of course, in the application to 
the Ising model we are discussing here one has n = 2, i.e. one works in the 
two-dimensional case; however the case of generic dimension n of the ambient 
space also makes sense from the mathematical point of view. The growth 
of a droplet and formation of new layers of the solid is described by the growth 
of profile u as m increases. 

The norm </'*(') here is related to the surface tension as follows. The surface 
tension 7^^(-) is a function defined over S^, the two-dimensional unit sphere 
in M^, and satisfying > and -/^^{-u) = -/^^{u) for all G §^ The 

surface energy of a closed surface C is defined to be 

H{M^)= [ ^^''{u,)ds, 
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where Ug is the unit normal to M"^ at s G M^. While defines the 3D shape 
of a crystal growing in space, the shape of monolayers growing on facets is given 
by the restricted 2D surface tension defined for n G by 



where the derivatives are taken at the "north pole" Ug = (0, 0, 1) G along 
all tangents G §^ to [9]. The function can then be extended to the 
whole by homogeneity of degree one, and </?*(•) is defined as the convex hull 
of thus defined 7^^(-)- However in the sequel (p* will be fixed, without any 
assumptions of smoothness or strict convexity (indeed one of the examples in 
Section [5] corresponds to a crystalline norm) . 

It is easy to see that the functional minimized in (jl.ip is the one-dimensional 
surface energy for the restricted surface tension. It turns out that minimization 
of this surface energy alone is sufficient to reconstruct most of the physics of 
monomolecular layers growth described in In particular, if ip*{-) is the 
Euclidean norm and a unit square, then as m grows, the first four monomo- 
lecular layers start as Wulff circles and then develop into "Wulff plaquettes" 
while from the fifth layer on all new layers appear as Wulff plaquettes identical 
to the underlying layers (Section [5]). 

In contrast, this simple variational model does not capture the thermody- 
namic fiuctuations, which render Wulff circles below a certain size unstable 
and prevent their formation for small m. Neither does it capture the micro- 
scopic (i.e., "finite- A^") structure of the Wulff plaquettes, whose boundaries 
are in fact separated with gaps that vanish in the continuous limit. A first- 
principle approach that takes proper account of these phenomena is due to 
R. Dobrushin, S. Shlosman and their coauthors and is presented in their works 



It is worth observing that a similar problem with additional restriction that 
u be a characteristic function of some set (i.e. that the droplet has exactly 
one layer) in the two-dimensional situation (i.e. when n = 2), the set is 
convex, and the norm 93* is Euclidean, has been studied in |16| . and for more 
general anisotropic norms (but for a somewhat different functional, namely, with 
penalization on the volume instead of the volume constraint), in [T3]. The latter 
problem will play an important role also in the present paper. Eventually, one 
has to mention that it is also very similar to the well-known Cheeger problem, 
the solutions of the latter being so-called Cheeger sets (see e.g. [H [121 [HI E] ) • 

Our aim in this paper is to study the variational problem (jl.ip in the two- 
dimensional case (i.e. when n = 2). This geometric optimization problem is 
considered without resort to the underlying lattice model or its continuous limit, 
allowing us to treat an arbitrary open domain Q and an arbitrary norm if* that 
is not necessarily strictly convex. In this setting we completely characterize 
the minimizers and the possible levels of u when the domain Q is convex. In 
particular it turns out that that except some degenerate situation, which can 
however happen only when is not strictly convex, the number of nonzero 
levels of u is at most two. 




[6l[T4l[71[8l[9l[T0]. 
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The basic tool we use is the auxiliary problem when u is a priori required 
to have a single nonzero level (i.e. is requested to be a characteristic function); 
namely, we show that in the two-dimensional case (n = 2) when Q is convex, 
the nonzero levels of solutions to the latter problem corresponding to different 
values of m as m grows can be arranged as a family of sets ordered by inclusion. 
Thus, solutions to problem (jl.ip can be seen as "towers" with levels solving the 
auxiliary problem. The assumption of convexity of is essential, as shown 
by a counterexample at the end of Subsection 14.11 The main result of the 
paper is formulated as Theorem 14.101 We conclude with an explicit example of 
solutions to (jl.ip for the case of a square Q = [0, 1]^ with the Euclidean norm 
and a crystalline norm. 

This work was inspired by some seminar talks of Senya Shlosman. After it was 
completed, we learned that a full description of the solutions to the variational 
problem (jl.ip when J7 is a square and (p* is generated by a physical Hamiltonian 
(in particular, is the Euclidean norm) has been independently obtained by him 
and loffe |15j by a rigorous continuous limit of a suitable lattice model. Their 
proof, together with an analysis of the microscopic structure of the solution 
and its behavior under thermal perturbations, will appear in the forthcoming 
publication [10]. 

2. Notation and preliminary results 

For a set E C M" we denote by \E\ its Lebesgue measure, by 1^; its charac- 
teristic function, by E its closure, by dE its topological boundary, and by E'^ 
its complement. 

In the following ip will denote the given (not necessarily Euclidean) norm 
over M". Given S C M" and x G M", we set 

distip{x, E) := inf ip{x — y), d^{x) := dist^p^x, E) — dist^p^x, E'^). 

The value (x) is the signed distance from x to dE and is positive outside E. 
Notice that at each point where d^ is differentiable one has (see [3j) 

(2.1) (/'*(V(i^) = 1, 1/ • Vd^ = 1 for all e dp*(yd^), 

where ip* denotes the dual norm of (p defined as 

■= max{^ • r/ : y?(r/) < 1} 

and dip* denotes the subdifFerential of ip* in the sense of convex analysis. In 
particular 

Vd^ 



where is the exterior Euclidean unit normal to dE. 

We define the anisotropic perimeter of a set E C as 
(2.2) 

P^{E) := supjy divT? dx: r] £ C7o^(M"), ipir]) < l| = y ^p* {i^^)dn''-\ 

where Cq(M") stands for the set of continuously differentiable functions with 
compact support is M", d*E is the reduced boundary of E according to De 
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Giorgi and T-O' stands for the /c-dimensional HausdorfF measure. We will usually 
identify a set E of finite perimeter with the set of its density points (i.e. points 
of density 1). 

Given an open set C M" we define the i?y-seminorm of f G BV{Q) as 



/ tf*{Dv) := sup \ V div 
Jn [Jn 



We let := {x\ <^(x) < 1}, usually called the Wulff shape, be the unit ball 
of If. Observe that -P(^(W(p) = n|W(^|. 

In the sequel, given x G M" and r > 0, we set Wr{x) := x + rW^p (a Wulff 
hall of radius r with center x). In this notation the reference to a norm is not 
retained for the sake of brevity, but always silently assumed. When 93 is the 
Euclidean norm, we will use a more common notation Br{x) instead of Wr{x) 
and P instead of P^ . 

Definition 2.1. Given an r > 0, we say that E satisfies the rW^- condition, if 
for every x G dE there exists an y & such that 

Wr{y) C E and x G dWr{y). 

Observe that, if E is convex, then E^ satisfies the rVF^-condition for all r > 0. 
We conclude the section by recalling the following isoperimetric inequal- 
ity HZ]. 

Proposition 2.2. For all E CW^ such that \E\ < +00 there holds 

n-l 

(2.3) PAE)>^^^PAW^)- 

\w^\ — 

3. Existence of minimizers 

Notice that, since the total variation is lower semicontinuous and the con- 
straints are closed under weak BV convergence, by direct method of the calculus 
of variations one immediately gets existence of minimizers of (jl.ip . 

Proposition 3.1. For any m > there exists a (possibly nonunique) minimizer 

of dni). 

For every u G L^{W^) and j G N we set 

(3.1) Ej := {u > j}. 

It is worth observing that whenever u{-) takes values in N, one has 

00 

(3.2) u = Y,^E, 

i=l 

and 

(3.3) / ip*{Dn) = J2PAE^)■ 



V dx = 1 
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Remark 3.2. It is worth observing that, if we let be a minimizer of (jl.ip 
for a given m > 0, then the normahzed functions Vm ■= Um/fn converge, as 
m — )• cxD, up to a subsequence, to a minimizer of the problem 

mini / (p*{Dv): V e BV(E."'),v = on n"", [ 

which is closely related to the so-called Cheeger problem in |12j . 

The following assertions hold true. 
Proposition 3.3. If u is a minimizer of (jl.ip . then u G L°°(M"). 
Proof. Assume by contradiction that \Ej\ > for all j £ N. Notice that 

lim \Ej\ = 

(since otherwise u would not be integrable). Given xq € Q we let 

Uj := mm{u,j) + 1wr^{xo) 
where the radius Rj is such that 

/ Uj = / u = m, 
Jn Jn 

that is (keeping in mind (j3.2p ) 

\W^\R^ = Y,\E^\, 

i>j 

and choose j € N big enough so that Wr^{xq) C ^l. 
Letting f{t) := n| so that 

PUWR^{xo)) = f{\W^\R]), 

we have 



n 



*{Duj) < [ ^*{Dmm{u,j))+P^{WR^{xo)) 
Jn 

= J ^*{Dmm{u,j)) + f{\W^\Rj) 

< / If* {D m.m{u, j)) + y ^ f{\Ei\) by the concavity of / 

< [ ip*{Dmm{uJ))+y2P^{Ei) by (USD 

= [ ^*{Du) by (13:31). 

Jn 

the second inequality being strict unless \Ei\ = for all i > j, > j, thus 
leading to a contradiction. □ 

Proposition 3.4. Lei C 6e star-shaped. Then Problem (jl.ip is equivalent 

to the following relaxed problem 

(3.4) 



mm 



/ Lp*{Du) : u G BV{W), u = on 0^ u{-) G N, udx>m 
Jn Jn 
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Namely, the minimum values and the minimizers are the same for both prob- 
lems. 

Proof. It is enough to show that any minimizer u of (|3.4p satisfies 

(3.5) / udx = m. 

Jn 

To this aim let be star-shaped with respect to xq and assume by contradiction 
that (|3.5p is violated. Let ux{x) := u{xo + A(x — xq)) for any A > 0, so that 
u\ G BV{MP), u\{-) G N, while, by star-shapedness of 17, one has u\ = outside 
of $7 for every A > 1. Then there exists a A > 1 such that ()3.5p holds with u 
replaced by ux- However 

/ ^*{Dux) = X^-'"" [ ip*{Du) = X^-"" [ if*{Du)< [ ip*{Du) 

Jn Jxo+x(n-xo) Jn Jn 

(the second equality is due to the fact that Q C xo + A(0 — xq) for A > 1, while 
u = outside of O), contradicting the minimality of u. □ 

4. The convex two-dimensional case 

In this section we shall assume that n = 2 and 17 C is a convex open set. 
Given E CM."^ and an r > 0, we define the set E"^ C E hy the formula 

f U {Wr{x) : Wr{x) CE}, if r > 0, 
(4.1) E^ := <^ 

[ E, if r = 0. 

Notice that, if ii^ is a convex set, then E^' is convex and satisfies the rWip- 
condition. The set E"^ is called the Wulff plaquette of radius r relative to E. 
The following assertion holds. 

Lemma 4.1. Let E C M? be a convex open set satisfying the rW^p- condition 
for some r > 0, then E = E'' . 

Proof. One has E^ C E. On the other hand, dE C dE^ because E satisfies the 
rWi^-condition. Minding that E, and hence is convex, we get E = E"^ . □ 

It is worth mentioning that convexity of the set E is essential in the above 
Lemma 14. li In fact, if A, B and C are the vertices of an equilateral triangle 
l\ABC with sidelength 1, then letting 

E := By^iA) U By^iB) U By^iC) U AABC 

we have that E satisfies the |VF<^-condition with respect to the Euclidean norm, 
but 

E'/'' = B,/^{A) U i?i/2(i?) U B,/,{C) + E. 
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4.1. Isoperimetric sets. We consider the constrained isoperimetric problem 
(4.2) min{P^(^) : E C n, \E\ = m £ [0, 

which corresponds to Problem (jl.ip under the additional constraint that u is a 
characteristic function. Clearly, the minimizers to this problem exist and the 
assertion of Proposition 13.41 remains valid for this problem. 

Let > be the maximal radius R such that Wr{x) C Q for some x £ 
Q, and let G [0, i?^] be the maximal radius r such that Q satisfies the 
rWi^-condition (we set for convenience := if does not satisfy any rW^- 
condition). Observe that in the Euclidean case one has 

1 

Tq = - — 

ll^llL°°(9f2) 

where k stands for the curvature of d^l. 

Lemma 4.2. Let m G (0, and let E he a minimizer of ()4.2p . Then E 
is convex and there exists an r > (depending on m) such that E satisfies 
the rW^-condition and each connected component of dE D is contained in 
dWr{x), for some Wulff hall Wr{x) C Vt, (with x depending on the connected 
component of dE n Oj. 

Remark 4.3. Recall that here and in the sequel when speaking of the properties 
of a set E of finite perimeter we actually refer to the respective properties of the 
set of its density points. In particular, a minimizer E of ()4.2p is not necessarily 
convex, but the set of its density points is (and hence, in particular, the closure 
E is convex). 

Proof. We divide the proof into four steps. 

Step 1. Let us first show the convexity of E. As in [U Theorem 2] we can 
uniquely decompose E as a union of (measure theoretic) connected components 
{-Ejjjg/, where / is finite or countable, such that 

1^1 = and P^{E) = Y,PAEi)- 

As in [2, Proposition 6.12], one shows by the isoperimetric inequality and 
the minimality of E, that the number of connected components is finite and 
the boundary of each connected component Ei is parameterized by a finite 
number of pairwise disjoint Jordan curves. In particular, the boundaries of 
two different connected components do not intersect. Further, using lemma 6.9 
from [2] one has that the perimeter P^{Ei) of a measure theoretic connected 
component Ei having the boundary parameterized by Jordan curves {Oj}^!^ 
(all parameterized, say, over [0, 1]) is given by 



T~, Jo 



where 'i/^ : — )• M is some convex and 1-homogeneous functions (in fact, ^p := 
if* oR^ R being the clockwise rotation of by vr/2, see corollary 6.10 from |2j). 
Hence, using Jensen inequality one shows that the convex envelope of Ei has 
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lower (anisotropic) perimeter than Ei itself, and minding that it also has greater 
volume (as well as the fact that the assertion of Proposition 13.41 is valid for 
Problem (j4.2p ). one has that each Ei is convex. 

Finally, if E is not connected, recalling that il. is convex we can translate 
a connected component inside 0, in such a way that its boundary touches the 
boundary of another connected component (this does not change neither the 
perimeter nor the volume), and taking the convex envelope of the resulting 
set we obtain again a set with greater volume and strictly lower anisotropic 
perimeter, hence a contradiction which shows that E is convex. 

Step 2. Reasoning as in |13^ theorem 4.5], where the authors consider the 
related problem mm{Py,{E) — X\E\ : E cil, X > 0} instead of ()4.2p . one gets 
that each connected component of dE D 0, is contained in dWr{x), for some 
X G and r > 0. 

Moreover, as in [2, theorem 6.19] one can show the existence of a (possibly 
nonunique) Lipschitz continuous vector field n: dE — such that n{x) G 
dip*{i'{x)) for "H^-a.e. x G dE. In particular div,-?! G L°°{dE), where divi-n := 
dr{n ■ r) denotes the tangential divergence of n (here and below r and u denote 
the Euclidean unit tangent and exterior normal vectors to dE respectively) and 
corresponds to the anisotropic curvature of dE (cfr. pTl [3]). 

Without loss of generality we may assume that divi-n is constant along every 
maximal segment contained in dE (if not, we can substitute n over the segment 
by a convex combination of its values on the endpoints of the segment; one 
would then still have n G d(p*{u) along the segment because u is constant there 
and dip*{-) is convex). In particular, if a connected component S of dE n O is 
contained in dWr{x), then n{y) = {y — x)/{rip{y — x)) for T-L^-a.e. y G S. 

Step 3. We now prove that E satisfies the rVF^^-condition for some r > 0. 
Since E is convex, it is enough to show that 



This follows by a local variation argument as in the proof of Lemma [4.91 b elow . 
Let us fix xi G S, where S is a connected component of dECii}, and X2 G dE\Ti. 
We know from the previous step that S is contained in dWr{x) for some x £M? 
and r > 0. We distinguish four cases. 

Case 1. There are two disjoint open sets C/j, i = 1,2, such that Xi G Ui 
and Ui n dE do not contain segments. Let ipi,ip2 be two nonnegative smooth 
functions, with support on Ui,U2 respectively, such that 



where i' stands for the exterior Euclidean unit normal to dE. We consider a 
family of diffeomorphisms such that 

^'(e, x) := X + eipi{x)n{x) — eip2{x)n{x) + o(e) 

for e > small enough. By (j4.4p . the term o(e) can be chosen in such a way 
that 

(4.5) \E'^\ = \E\ for all e > small enough. 



(4.3) 



div,-?! < — 
r 



n^-a.e. on dE. 



(4.4) 
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with := ^{e,E) C il. We then have 

P^{E') = P^{E) + - f il,i{z)^*{u{z))(m}{z) 
JUindE 

-el 'ilj2{z)dwrn{z)^*{v{z))(m}{z) + o{£). 
Ju2r\dE 

As e — 7- O"*", by minimahty of E, we get 

-/ Mz)v*Hz))dn\z)> [ Mz)divrn{z)ip*{u{z))dn\z) 
JUindE Ju2ndE 

which in view of (j4.4p gives (j4.3p . 

Case 2. We can find two maximal segments ii,£2 C dE such that Xi G ii, and 
we define E"^ by shifting ii by cie parahel to itself outside E, and by shifting 
£2 by C2£ inside of E, with ci, C2 so that ()4.5p holds, that is 

(4.6) Ci\ii\ = C2\e2\. 
By |13t Lemma 4.4] we have 

P^{E^) = Pip{E) + ciaie - 02026 + o(e). 

where ai, Q2 are respectively the (Euclidean) length of the face of Wy, par- 
allel to ii, £2- By minimality of E, letting e ^ we obtain cioi > 0202- 
Recalling (|4.6I1 . we finally get 

- = 1^ > 1^ = div^n(z) for z G £2- 
r \£i\ K2I 

Case 3. There is a maximal segment £1 C and an open set U2 such that 
xi £ £1, X2 £ U2 and C/2 H does not contain segments. We proceed by 
combining the previous strategies and we define the set E^ by shifting £1 by 
£ parallel to itself outside E, and then taking the image of the resulting set 
through the diffeomorphism 

^'(e, x) := X — eip2{x)n{x) + o(e), 

where ip2 is a nonnegative smooth function supported on U2 satisfying 

(4.7) / M^)^*H^))dnHz) = \£,\. 
JU2ndE 

This condition guarantees that the volume change after these two operations 
is of order o(e), so that the extra term o(e) in the definition of ^ is chosen in 
such a way that ()4.5p holds. Reasoning as above, we get 

P^iE') = P^{E) +aie-e [ Mz)divrn{z)ip*ii^{z)) dn^z) + o(e), 

Ju2ndE 

which gives, by minimality of E, 

ai = ^> / Mz)diVrn{z)ip*{u{z))d'H\z) 
Ju2r\dE 

which gives (j4.3p . recalling (j4.7p . 

Case 4. There is a maximal segment £2 C dE and an open set Ui such that 
xi G Ui, X2 G £2 and Ui Ci dE does not contain segments. This case can be 
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dealt with reasoning as in the previous case, by shifting £2 by e inside E and 
defining 

^'(e, x) := X + eilJi{x)n2{x) + o(e). 

Step 4. From (14. 3p it foUows that the radius r in Step 3 does not depend 
on the connected component S. In particular, every connected component of 
dE n is contained in dWr (x) , for a fixed r > (while x depends in general 
on the connected component). □ 

Consider now the function v{r) := It is clearly constantly equal to 

\Q\ for r < and to zero for r > Rq, while over [r^,RQ] it is continuous 
and monotone decreasing. In particular, for all m G there exists a 

unique value G [rn, Rn] such that v{rm) = m. 

From the isoperimetric inequality (j2.3p and Lemma [4. 21 we get the following 
statement. 

Proposition 4.4. Let C 6e convex, and let E be a minimizer of (14. 2p 
with m £ [0, Then 

(a) either E = J^^"", if m > \Q^^\, 

(b) or E is the closure of some convex union of Wulff halls of radius Rq, if 
mG_[i?2|Ty^|,|0««|], 

(c) or E = W y ^yi^j/ for some x £ ^l, if m < i?^|W^|. 

Proof. We can assume m G (0, By Lemma 14.21 there exists an r > 

(depending on m) such that E is the closure of a union of Wulff balls of radius 
r, hence E aVf and r < 

If m > then necessarily r < R^ and E = Q^', since otherwise we could 

find a connected component of dE n Q which is not contained in the boundary 
of a Wulff ball, contradicting Lemma 14.21 In particular, we have r = r^. 

If m G [i?^|VF^|, then r = Rq, since otherwise E would coincide with 

the set Cl^ (with r < Rq) which has volume strictly greater than 

If m < the result follows by the isoperimetric inequality (j2.3p . □ 

Remark 4.5. It is worth noticing that, if is strictly convex, then there exists 
a unique Wulff ball Wr^{x) C O, and thus = Wji^{x). In other words, 
the case (b) of the above Proposition 14.41 reduces to case (c). Therefore, either 
E = Cl''^, if m > OT E = W^ ^i^^ ^ [x) for some x G if m < \VL^^\. 

We now state an easy consequence of Proposition 14.41 showing that solutions 
to problems (j4.2p with decreasing volumes may be arranged as a decreasing 
sequence of sets. 

Corollary 4.6. Let he convex and let rrij he a decreasing sequence such that 
rrij G (0, for all j. Then, there exists a sequence of sets Ej such that 
Ej^i C Ej C \Ej\ = rrij and each Ej is a minimizer of (14.21) with m := rrij. 

Note that the convexity assumption of the set O is essential in the above 
result. In fact, reasoning as in [ TTl section 6] with the example of $7 a couple of 
circles connected by a thin tube (like a barbell considered in |1H section 6]), one 
provides a family of minimizers of (14. 2p with decreasing volumes which cannot 
be arranged as a decreasing sequence of sets (see Figure 14.11 below) . 
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Figure 1. C nonconvex (two circles connected with a thin 
tube) and two minimizers of (|4.2p which cannot be included one 
into another. 



4.2. Isoperimetric towers. We return now to the original problem (jl.ip . 
Here and below we let u £ L^(]R^) be an arbitrary minimizer of this prob- 
lem and Ej be its level set corresponding to a j G N, as defined by (|3.ip . The 
following result follows directly from Corollary 14.61 

Proposition 4.7. If Q is convex, then for all j G N the set Ej is a minimizer 
of problem (j4.2p with m := \Ej\ (in particular Ej is convex). 

Proof. If the assertion is not true, then considering a sequence of sets Ej of 
minimizers of (|^ (with m := \Ej\) such that Ej_^_^ C Ej C ri, \E'-\ := \Ej\ 
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(the existence of such a sequence is guaranteed by Corollary I4.6P , and setting 

j 

we get 

the strict inequality being due to the fact that one of Ej is not a minimizer 
of ()4.2p (with m := \Ej\) by assumption. On the other hand, 

u' dx = I udx = m, 
Jn 

since the level sets of u' and u have the same volume by construction. This 
would mean that u is not a solution to problem (jl.ip . □ 

Remark 4.8. Observe that, by Proposition 14.71 and Lemma 14.21 each set Ei is 
convex and each connected component of dEi n ^2 is contained in dWr^{xi) for 
some Wulff bah Wr^{xi) C il. 

Lemma 4.9. Let Si, Sj be connected components of dEi n Q and dEj n O 

respectively, with j > i, such that 

c dWr^ix^) c n 
Sj c dWrjixj) c n 

(4.8) ^(Si-Xi)c -{Sj-Xj) 

n rj 

for some Xi,Xj £ Q, ri,rj > 0. Then ri > rj. 

Proof. It is enough to consider the case j = i+1. We can also assume Si ^ 5'j+i, 
otherwise there is nothing to prove. As in Figure 14.21 there are two cases to 
consider. 

Case 1. There are two points yi £ Si, y^+i € 5j+i and two open sets Ui C 
and f/j+i C ft such that yi £ Ui, yi+i G C/j+i, Ui n C/j+i = and Ui n Si as well 
as C/j+i n S'j+i does not contain segments. Consider a smooth function ipi with 
support on Ui. It generates a one-parameter family of diffeomorphisms of Ei 
defined by 

^i{e, x) := X - e^i{x)ni{x) 
for all sufficiently small e > 0, where 

ni[x) := -. 

riip[x - Xi) 

Consider now a one-parameter family {^'j4.i(e, •)} of diffeomorphisms of i^j-fi 
such that ^'j+i(0,x) = x for all x G i^j+i, ^'i+i(e, •) — Id is supported in C/j+i 
for all e > 0, while 

^fj+i {e,x) := x + eipi+i (x)nj+i (x) + o{e) 
as e — )• 0^, where ipi+i is some smooth function (with support in f/j+i), and 

X - Xi+l 



ni+i{x) :-- 



ri+iip{x - Xi+i) ' 
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We choose ^i+i so that the sets Ef := ^i{e,Ei) and E?_^_-^ := \I'j+i(e, 
satisfy 

\Et\ + \E!^,\ = \E,\ + \E,+,\ 

for ah sufficiently smah e > 0. Denote by I'j the exterior Euchdean unit normal 
to dEj. Since 



\Ef\ = \Ei\-e [ iP,{z)^*{iy,{z))dn\z)+o{e), 
JdEiDUi 

lEf+il = \Ei+^\ + e I i^i+i{z)ip*{vi+i{z)) dn\z) + o{e), 



as e — )• 0^, we have 

^|;,{z)cp*ily^{z))dn\z) 



(4.9) 



Letting now 



dE.nUi 

ipi+i{z)ip* (i^i+iiz)) dn^iz). 

dEi+inUi+i 



Ue:=U- 1e, - + lEf + lEf^^ = J] ^^fe + ^^f + 

we have /^^ dx = J^u dx for all sufficiently small e > 0. Recall that 



-^,{z)^*{u,{z))dn\z) 

+ e [ -^iJi+i{z)ip*{i^i+i{z))dn\z) 



+ o(e). 

As e — )• 0"*", by minimality of li, we get 

-- [ M^)ip*{iyi{z))dn\z 

JdEiHU^ 

1 



+ / ^P^+l{z)ip*{lJi+^{z)) dn\z) > 0, 

JdE,+inu,+i 



which together with (j4.9p implies the thesis. 

Case 2. We can find two maximal line segments ii C Si and C -Sj+i- We 
define then E? by shifting the segment ii by Cjff parallel to itself inside Ei and 
-E-f+i by shifting the segment ^j+i parallel to itself outside of by Q+ie with 
Cj and so as to satisfy 

\Ef\ + \Ef^,\ = \E,\ + \E,+^\ 

for all e > sufficiently small. Since 

\Ef\ = \Ei\-Ci\i^\e + o{e), 

= \Ei+i\ + Cj+i|£i+i|e + o(e), 
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as e — )■ 0"^, we have 

(4.10) c^\ii\ = Ci+i\li+i\. 

Letting again, as in Case 1, 

Ue:=u- 1e, - 1e,+, + Iei + = Yl + Ie; + Ie^^, , 

we have Ug dx = f^u dx for all sufficiently small e > 0. On the other hand, 
by \13\ lemma 4.4] 

/ ip*{Due)= / ip*{Du) - aaiE + Ci+iai+ie + o{e). 
Jn Jn 

where a^, a^+i are respectively the (Euclidean) length of the face of parallel 
to ii, ii+i- By minimality of u, letting e — )• 0'*' we obtain CjOj < Q+iaj+i. 
Recalling (|4.10p . we get = \ii\/ai > \£i+i\/ai+i = rj+i. 

Notice that in this proof we do not have to deal with the situation depicted 
in Cases 3 and 4 of the proof of Lemma 14.21 due to condition (j4.8p . In fact, 
the latter implies that if Si contains a line segment ii, then the line segment 
ij := Xj + {li — Xi)rj/ri is contained in Sj. Otherwise, if there is a neighborhood 
Ui of a point of Si such that Si D Ui does not contain any line segment, then 
Uj := Xj + {Ui — Xi)rj/ri is a neighborhood of a point in Sj such that Sj (1 Uj 
does not contain any line segment. □ 



We are now able to prove the following result giving the complete character- 
ization of solutions to problem (jl.ip . 

Theorem 4.10. Let 17 C be convex and set jj^a^^ := ||m||oo- Then one of the 
following cases holds. 

a) There exists an f ^ [tq, Rq) such that Ej = for all j < jmax- IiT' this 
case 

U = Jmaxljl'' 

(in particular, if f = r^, then u = jma.x'^n)- 

b) There exists an f £ {rQ,Rn) such that Ej^^^ = Wf(x) for some x G O 
with Wr{x) C , and Ej = Cl^ for all j < jmax- Ii^- this case 

U = ^Wf{x) + (imax - l)lf7*'- 

c) There exists an f £ {0,r^] such that Ej^^^ = Wf{x) for some x G Q 
with Wr{x) C il, and Ej = Q for all j < jmax- In this case 

U = 'i-Wr(x) + (jmax - l)ln 

(note that this condition may hold only when tq > 0). 

d) Every Ej is the closure of a convex union of Wulff halls of radius Rq, 
for all j < j 

max- 

Remark 4.11. Observe that case d) of Theorem 14. 101 is the only case where the 
number of nonzero level sets of the minimizer may be bigger than two. 
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Proof. We may assume jmax > since otherwise the result fohows directly 
from Proposition 14. 4i 

By Remark 14.81 for all i < jmax the set Ei is convex and each connected 
component of dEiCiCl is contained, up to a translation, in dWr^ (xi) for some rj > 
0, Xi G fl. Moreover, if dEi n Q and dEi^i n Q are nonempty, from the inclusion 
Ei^i C Ei it follows that we can always find two connected components Si C 
dEi n Q and /Sj+i C dEi^i n il. satisfying the assumptions of Lemma 14.91 By 
Lemma urn] we then get rj > rj+i for all i < jmax- 

Recalling Propositions 14.71 and 14.41 this leaves only the following possibilities: 

(i) Ei = Cl^'- and -Ej+i = fi^'+i, with rj > rj+i. In this case, if rj > tq so 
that Cl^''- 7^ Cl, then = r^+i, hence Ei = Ei^i = while if r^+i < 
< fn we have Ei = Ei^i = Cl, and we may just set = rj+i := so 
that we still have Ei = Ei^i = Cl^'^ = Q. 

(u) Ei = and Ei+i = Wr.+Jxj+i), with > r^+i. 

(in) £'i = Wriixi) and E'j+i = ly^.+i (a^i+i), with > r^+i. 

(iv) Ei is a closure of a convex union of Wulff shapes of radius Rn and 
Ei^i = (xj+i), with Rq > Tj+i (note that the case where Ei^i 

is the closure of a convex union of Wulff balls of radius Rq and Ei = 
Wri{xi) with Rfi < ri is impossible). 
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(v) Ei and -Ej+i are both the closure of a convex union of Wulff shapes of 
radius Rq. 

Thus there is J G {0, . . . , jmax} and f G [tqjRq) such that, for every i < j, 

(A) either Ei = , and Ei = Wri{xi) with ri < f for all i > j; 

(B) or £'j is the closure of a convex union of Wulff shapes of radius R^, and 
Ei = Wri{xi) with Tj < i?f7 for all i > j. 

We want to show that j = jmax — 1 or j = jmax- 

Assume by contradiction that jmax > J + 2. As rj+i < f in Case A, and 
Tj < R^ in Case B, without loss of generality we can assume that Ej^i = 
W^T-j+i is contained in the interior of Ej. Then, we can choose e,e' > 

sufficiently small so that, letting r'j^^^^ := rj^^^ — e' and rj_^-^ := r^+i + e, we 
have Wr^^^{xj+i) C Ej^^^, and 

(4-11) ^'L. + ^1+i = d.. + ^|+i- 

From (I4.1ip we then get 

e' = ^e + o{e), 

Jmax 

and hence 

J+1 I Jmax Jmax 

= l-e ^^-+^~^^— 

where the error term o(e) is negative whenever r^+i = rj^^^. Representing u as 
u = u + 1 + Ik , and letting 

J+1 Jmax 

we then get 

! Jn 

for sufficiently small e > 0. Since f^2udx = f^2u'^dx, this contradicts the 
minimality of u. One has therefore that j = jmax — 1 or j = jmax, which 
concludes the proof. □ 

5. An EXPLICIT EXAMPLE 

5.1. A square with the Euclidean norm. Let now Q := [0, 1]^ and let if 

be the Euclidean norm on M^. From Theorem 14.101 we obtain the following 
characterization for the minimizers of (jl.ip . 

Proposition 5.1. Let Q and ip he as above. 

(i) If m £ {n — l,n7r/4), with 1 < n < A, we have jmax = n, Ej^^^^^ = 
Br{xo) C and Ej = Q,^ for j < jmax, with 

r = \J (n — m— l)/(4(n — 1) — nvr) 

(ii) If m €z [mT/4,n], with 1 < n < 4, we have jmax = n cL^d Ej = for 
j < jmax, with r = ^{n- m)/((4 - 7r)n). 
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(iii) If m> A we have 



(5.1) 



Jmax S 



2 + 
2^ 



-m 



2 + 
2^^ 



-m 



+ 1 



and £'j = n'' for j < jmax, with r = {1 - m/j^^)/{A - tt). 

Proof. Clearly, = 0, Rn = 1/2. By Theorem 14.101 for all m > we have one 
of the following two possibilities. 

Case A. Ej = for all j < jmax with 

m = jmaxlf^n = Jmax (l " (4 - ^)r2) r G [0, 1/2]. 
It then follows 



r = r^Omax) : = 
and EfcT PiEj) = FA(imax), where 



m 



(4 - vr)i. 



Notice that 



F'^{x) =4-^/4^- 



2x — m 



^/x{x — m) 

which implies that Fji{x) is increasing for x > ^^'^^ m and decreasing for m < 



X < ^j^-iTT'- As a consequence we have 



jmax G {j^,J^ + 1}, where := 



2 + 
2V^ 



Case B. -B^^^,^ = Br{xQ) C and Ej = O*" for ah j < jmax with 
m = vrr^ + (j^ax - ^)\^''\ 



1) 1 



Jmax \ 2 
TT r 



re (0,1/2). 



It follows that 



Jii 



1 — m 



(4 - 7r)(jmax - 1) - vr 



and Ej=T ^(^j) = ^B(jmax), where 

Fb{x) := (x-l)P(!^'^s(^)) + 27rrB(x) 



= 4(x - 1) - 2V4 - vr , 
Notice that the derivative 



4 - vr 



(x — 1 — m). 



F^(a;) = 4-^4^- 



2(x — 1) — m 



TT 



4 - TT 



4 - TT 



(x — 1 — m) 
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Assuming x > jr^, we then have that Fb is increasing for 

2 + /^ 4 + 0F 
X > „ _ m + 



2^ 2(2 + V^) 

and decreasing otherwise, so that 



jmax G + where 



■m + 



2^ 2(2 + V^) 



as soon as jmax ^ 5. 

Observe that, if vr < m < 57r/4, then Case B cannot occur. In fact, supposing 
the contrary, we would have jmax G {j^ + 1} = {4, 5}, while one should have 

4 - 7rjmax/(jmax " 1) > 0, i.e. jmax > [4/(4 - 7r)J = 5 which implies jmax = 5. 

On the other hand either 

(1) jmax - 1 > vr/(4 - tt) with jmax >m+l, which would mean jmax > 6, 
so that this case is impossible, or 

(2) jmax - 1 < 7r/(4 - vr) with jmax <m+l, which means jmax = 4. 
This gives the contradiction proving that Case B does not occur. 

If m > 57r/4, we have to determine which one between Case A and B 
is energetically more convenient. By a tedious calculation one proves that 
min{ Fb ( j^ ), FB(j'^ + 1)} > min{FA(j^), ^^(j^ + 1)}, which shows that Case 
B still cannot occur. This shows claim (iii) of the statement being proven. 

Finally, with m < vr one can only have jmax £ {15 2,3,4}. It is easy to 
verify then that when m £ (n — l,n7r/4) and j = n G {1,2,3,4}, one has 
PbU) < FaU), so that Case B occurs, thus proving claim (i), while when 
m G [n7r/4, n] and j = n & {1,2,3,4}, one has Fa^J) < Fsij), and hence 
Case A occurs, thus proving claim (ii). □ 

Notice that is the volume of the (unique) Cheeger set Cq of Q, so that 
Proposition 15.11 implies that the functions Um/jmax converge to the character- 
istic function of Cq, according to the Remark [ 



5.2. A square with a crystalline norm. Now we set = [0, 1]^ as above 
and fii') = 1 1^1 1 + |f2|. Notice that <^ is a crystalline norm with Wulff shape 
= {{x,y) £ : |x| + \y\ < 1}. As before, we are able to characterize 
completely the minimizers of (II. ip . 

Proposition 5.2. Let Q and ip he as above. 

(i) If m £ (0,1/2], we have jmax = 1 CLnd Ei = Wrixo) C ^, with r = 

(ii) If m £ [1/2, 1), we have jmax = 1 and Ei = 0'', with r = y^(l — m)/2. 

(iii) If m = 1, then either jmax = 1 o^nd Ei = Q., or jmax = 2, Ei = fi'" and 
E2 = Wr{xo) C with r £ (0, 1/2]. 

(iv) If m> 1, we have 



(5.2) jmax G 



l + \/2 
m 



l + ^/2 
m 



+ 1 



and Ej = n'' for j < jmax, with r = y^(1 - m/jmax)/2. 
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Proof. The proof is similar to the one of Proposition 15.11 

Clearly, if m < 1/2, then jmax = 1 and Ei = Wr{xo) C 0, since the (rescaled) 
Wulff shape solves the isoperimetric problem. By Theorem l4.101 for all m > 1/2 
we have one of the following two possibilities. 

Case A. Ej = Q^' for all j < jmax with 

m = jmaxlJ^n = jmax (l " 2r^) T G [0, 1/2], 

which gives 

1 m 



x"^ — mx 



Jma 

and J2f:.r PA^j) = ^^(imax), where 

= x(4 - 4r) = 4x - 4 ^ 

Since the function Fa is increasing for x > ^~^^ m and decreasing for m < x < 
^^^ m, we have 



G + where / := 



1 + V2 
rn 



Case B. Ej^^^ = Wr{xQ) C and Ej = W'' for ah j < jmax, with r £ (0, 1/2] 
and 

m = 2r2 + (jmax - 1) (1 - 2r2) , 
and hence m > 1 because jmax 

> 2 and r < 1/2. 
If m = 1 then jmax = 2 and we can take any r G (0, 1/2]. 
If m > 1 then jmax > m + 1 and we get 



' jmax 1 



y 2 (jmax 2) 

and Ej=T PA^j) = FBijmax), where 

Fb{x) = 4{x - 1) -4{x - 2). 



I x — m — 1 
2{x - 2) 



Since the function Fb is increasing for x > ^^^ m + ^ and decreasing 
otherwise, we have 



e {j^,j^ + l}, where j^ 



l + \/2 3-^/2 
m H 



As in the proof of Proposition 15.11 when m > 1 we have to determine 
which one between Cases A and B is energetically more convenient. Since 
min{FB(j^), FbU^ + 1)} > min{FA(j^), FA{j^ + 1)} (again, by a tedious cal- 
culation like in the example with Euclidean norm), it follows that Case B can 
never occur. □ 
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